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Quantum Monte Carlo and density-matrix renormalization group methods are used to study the 
coupled spin-pseudospin Hamiltonian in one-dimension (ID) that models the charge-ordering insta- 
bility of the anisotropic Hubbard ladder at quarter filling. We calculate the temperature dependence 
of the uniform spin susceptibility and specific heat as well as the spin and charge excitation spec- 
tra of the system. We thereby show that there is a parameter and temperature region where the 
spin degrees of freedom are separated from the charge degrees of freedom and behave like a ID 
antiferromagnetic Heisenberg model, and that, outside this parameter region and above a crossover 
temperature, the spin excitations are largely affected by the charge fluctuations. We argue that 
observed anomalous spin dynamics in the disorder phase of a typical charge-ordered material a'- 
NaV2 05 may possibly be a consequence of this type of spin-charge coupling. 
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I. INTRODUCTION 

Charge-ordering (CO) instability has recently been one 
of the major topics in the field of strongly correlated elec- 
tron systems. Here, elucidation of the observed anoma- 
lous behaviors of electrons associated with the CO phase 
transition has been the central issue. The issue includes 
questions on the slow charge dynamics above the transi- 
tion temperature Tco as well as on the CO spatial pat- 
terns realized below Tco- A well-known example is the 
vanadate bronze a'-NaX^Os where the system may be 
modeled as a lattice of coupled ladders (or a trellis lat- 
tice) at quarter filling. 1 ' 2 ' 3 ' 4,5 Strong intersite Coulomb 
interactions between electrons are believed to be the ori- 
gin of the CO instability,^ In this material, the CO with 
a zigzag ordering pattern is observed below Tco = 34 
K)SiL2i2ii& and associated with this, a number of anoma- 
lous behaviors, which can be related to the slow dy- 
namics of charge carriers (or charge fluctuations), have 
been observed above Tco i 10 ! 11 : 12 ! 13 : 14 ! 15 : 16 ! 17 Anomalous 
response of the spin degrees of freedom has also been 
noticed iSiiSil^fliSi It seems therefore natural to wonder 
how, in such systems, the spin degrees of freedom be- 
have near the CO phase transition when they are on the 
slowly fluctuating charge carriers. In this paper, we con- 
sider this issue: i.e., what are the consequences of the 
slow charge fluctuation at T > Tco to the behavior of 
the spin degrees of freedom of the system? 

One of the simplest models that allow for such situa- 
tion is the anisotropic Hubbard ladders at quarter filling 
with the strong intersite Coulomb repulsions. We here 
use an effective Hamiltonian written in terms of the spin 
and pseudospin operators 5 * 1 ^* 2 ^ 3 . (where the latter rep- 
resents the charge degrees of freedom). This Hamilto- 
nian is derived from the Hubbard ladder model by the 
perturbation theoryf 5 *^ 2 ^ where the hopping parameter 
between the rungs of the ladder is assumed to be small 
compared with the onsite and intersite Coulomb repul- 



sions as well as the hopping parameter of the rung (i.e., 
the anisotropic ladder) £ Although the long-range CO is 
not realized in this model at T > (since it is the ID 
quantum-spin model), we can simulate anomalous behav- 
iors of the spin degrees of freedom under the influence 
of strong charge fluctuations. We will apply the quan- 
tum Monte Carlo (QMC) method to this model to cal- 
culate the temperature dependence of the uniform spin 
susceptibility and the spin and charge excitation spectra, 
thereby clarifying consequences of the interplay between 
its spin and charge degrees of freedom. The density- 
matrix renormalization group (DMRG) method with the 
finite-temperature algorithm will also be used to calcu- 
late the temperature dependence of the specific heat of 
the model. 

In this paper, we will first confirm that, in this cou- 
pled spin-pseudospin model, the spin exchange interac- 
tion is necessarily associated with the charge excitation; 
i.e., the spin excitations cannot occur without making 
the exchange of the pseudospins. We will then show 
that, nevertheless, there is a parameter and temperature 
region where the spin degrees of freedom behave like a 
ID antiferromagnetic Heisenberg model; i.e., the spin de- 
grees of freedom are 'separated' from the charge degrees 
of freedom in this region. We will moreover show that 
the spin system behaves in different manner depending on 
whether the temperature T is below or above a crossover 
temperature T* that is related to the pseudospin exci- 
tations; at T < T* , it behaves like a ID antiferromag- 
netic Heisenberg model with a T-independent effective 
exchange coupling constant J e ff with the large renormal- 
ization, whereas at T > T*, J e g decreases rapidly with 
increasing T, where the effective Heisenberg-model de- 
scription ceases to be valid. Because the parameter val- 
ues for a'-NaV205 are outside the region where the spin- 
charge separation is complete, we will argue that some 
experimental data may possibly be interpreted as conse- 
quences of this type of spin-charge coupling. 
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This paper is organized as follows. In Sec. II, we de- 
fine the coupled spin-pseudospin model that describes the 
spin and charge degrees of freedom of the anisotropic 
Hubbard ladder at quarter filling. Some details of the 
method of calculation are also given. In Sec. Ill, we 
present results of calculation, including the staggered sus- 
ceptibility for pseudospins, the spin and pseudospin ex- 
citation spectra, and the temperature dependence of the 
uniform spin susceptibility and specific heat. Discussion 
on the experimental relevance to a'-NaV20s and sum- 
mary of the paper will be given in Sec. IV. 



II. MODEL AND METHOD 

Our effective spin-pseudospin Hamiltonian for the 
anisotropic Hubbard ladder at quarter filling may be 
written as a sum 

H = H q + Hst (1) 
of the quantum Ising Hamiltonian for pseudospins 

L L 

i=l i=l 

and the spin-pseudospin coupling term 
L 1 

n ST = J 2 J2 ( s * • s *+i - 4) ( T + T i+i + H - c -) ■ ( 3 ) 

i=l 

The standard notation is used here. Sj and Tj are, re- 
spectively, the spin and pseudospin operators of spin- 1/2 
at site i, where T* = —1/2 (+1/2) means the electron is 
on the left (right) site on the rung of the ladder. L is the 
system size and periodic boundary condition is assumed. 
J\ is the energy scale of the pseudospin system and J2 is 
the coupling strength between the spin and pseudospin 
systems. 

The effective Hamiltonian Eq. (1) may be obtained 
from the second-order perturbation theory we have 
the relations J\ = 2Vj| and J2 = 4t|/Vj_, where i|| and Vy 
(t± and Vj_) are the nearest- neighbor hopping parameter 
and Coulomb repulsion of the leg (rung) of the ladder, 
respectively. We should then have J\ > J2, which we 
assume throughout the present work. We also assume 
the onsite Coulomb repulsion to be U — > 00. Relative 
strength of the transverse field applied to the pseudospins 
is measured by g = it±/Ji = 2t±/V^ . Note that g in the 
quantum Ising model represents the relative strength of 
the fluctuation of a charge in the rung: if we assume 
one electron in a rung, we have the prefactor gJ\/2 in 
the first term of Eq. (2), which is the difference between 
the energies of the bonding and antibonding levels of the 
rung, 2t±. Thus, if g (or t±) is large the electron is stable 
in the bonding level of the rung, but if g (or t±) is small 
the effect of Vji easily leads the system to CO. 



We use the conventional world-line QMC method for 
the analysis of the model. We use a 32-site cluster (where 
a site contains a spin and a pseudospin) with periodic 
boundary condition; the cluster-size dependence of the 
calculated results are examined by using clusters of up 
to 96 sites but we find no significant size dependence in 
the results. Because the model does not conserve the to- 
tal pseudospin, we have examined a number of ways of 
the spin flips and confirmed that available analytical re- 
sults are reproduced correctly^* The maximum-entropy 
method is used to calculate the dynamical quantities 
like the spin and pseudospin excitation spectra. The 
DMRG method with finite-temperature algorithm'" is 
also used for the calculation of the temperature depen- 
dence of the specific heat of our model; the method en- 
ables us to access to sufficiently low temperatures. 

III. CALCULATED RESULTS 

A. Staggered susceptibility for pseudospins 

We first consider the nonlocal spin susceptibility de- 
fined as 

xa = Aa ((sj(-iA)sf) - (s;><sf» (4) 

Jo 

where Sj(—iX) is the Heisenberg representation of Sj 
and (• ■ ■) is the canonical average. Xij 1S Fourier trans- 
formed to the q-dependent susceptibility x(o)i which we 
calculate by the QMC method; the q — > limit gives 
the uniform spin susceptibility x(T) and the staggered 
spin susceptibility is defined as x(q) at <? = 7r. In the 
following, we calculate the susceptibilities for spins and 
pseudospins, whereby we use the subscripts S and T as in 
Xs(l) and Xt(q), which stand for the susceptibilities of 
the spin and pseudospin degrees of freedom, respectively. 

The phase diagram of the quantum Ising model TIq (Ti. 
at J2 = 0) is well known; 27 at T = there is a long-range 
order for g < 1 (g = 1 is a quantum critical point), which 
corresponds to the zigzag (or 'antiferromagnetic') CO. 
The calculated staggered susceptibility for pseudospins 
is shown in Fig. 1, where we find that it shows divergent 
behavior at T — > for g < 1. The dispersion relation 
of the pseudospin excitation observed in the calculated 
dynamical structure factor (shown in Fig. 2, see below) 
agrees well with the exact result^ 

uj q = -y \/l + 9 2 + Zgcosq. (5) 

We find in Fig. 1 that the inclusion of the coupling term 
Ti-ST, which introduces the quantum fluctuation via the 
factor T^TJ , suppresses the divergence. Thus, the in- 
clusion of the spin degrees of freedom via Hst tends to 
suppress the instability to the long-range order of pseu- 
dospins. 
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FIG. 1: Temperature dependence of the staggered suscepti- 
bility for pseudospins Xt(vt) calculated for the coupled spin- 
pseudospin Hamiltonian. 



B. Spin and pseudospin excitation spectra 

The dynamical pseudospin structure factor St(<7,w) is 
denned as 

ST(q,r) =±J2 e ~ iqirj ~ ri) {T r z i (T)T r * j (0)) (6) 



1 f°° 

St(q,t~) = — / duj ST(q,u))K(u),T) 
t Jo 



(7) 
(8) 



where St (<7, t) is the Fourier transform of the imaginary- 
time correlation function. We use the maximum entropy 
method for the inverse Laplace transformation (or an- 
alytical continuation) to obtain St(q, oS) from St{q,t). 
The dynamical spin structure factor Ss(q,oj) is similarly 
defined by replacing the pseudospin operator T* with the 
spin operator S*- 

The calculated results for the pseudospin excitation 
spectra at low temperature (k^T = O.U2) are shown in 
Fig. 2, where we find that the spectra are under strong 
influence of the spin-pseudospin coupling term J 2 . With 
increasing the coupling strength J2/J1, the peak of the 
pseudospin spectra shifts to higher energies and simul- 
taneously the spectra are broadened. Thus, the lower- 
energy edge of the peak is not affected strongly by the 
coupling strength J 2 , at least when g is large. We sup- 
pose that the scattering of the pseudospin excitations due 
to spin excitations causes the broadening of the spectra. 



The calculated results for the spin excitation spectra at 
low temperature are shown in Fig. 3, where we find that, 
in contrast to the pseudospin spectra, the spin excitation 
spectra change very little; i.e., the peak position, width, 
as well as the shape of the spectra are not affected by 
the parameter J\ when g > 1. When g is small, however, 
the peak position is slightly shifted to lower energies with 
increasing the value of Ji (see Fig. 3 (a)). 

The dispersion relation of the spin and pseudospin ex- 
citations calculated at low temperature are summarized 
in Fig. 4, which are obtained as the momentum depen- 
dence of the peak position of the spectra. For compari- 
son, we show the dispersion of the quantum Ising model 
in Figs. 4 (a) and (c); the gap opens when g > 1, which is 
closed at q — 7r when g — ► 1, leading to the 'antiferromag- 
netic' long-range order (or zigzag CO)i^^ We note that 
the gap remains open irrespective of the value of g when 
we include the coupling term J 2 . In the left panels and 
right panels of Fig. 4, we present the same dispersion re- 
lations uj q , but in a different energy scales, i.e., uj q / J\ and 
L/q/ J 2 . We find that, unless g is small, the spin excita- 
tion spectra are always inside the charge gap, i.e., inside 
the gap of the pseudospin excitation spectrum. Thus, 
when the charge gap is large, the energy scale of the spin 
excitations is separated from the high-energy charge ex- 
citations. With decreasing g, however, the energy of the 
charge excitation decreases at the momentum q = n to 
couple with the spin excitations. 

In Fig. 4 (b) and (d), we find that, for g > 1, the 
dispersion of the spin excitation spectra scales very well 
with J 2 ; i.e., it does not depend on the value of J\. More 
quantitatively, the dispersion of the calculated spin exci- 
tation spectra is fitted well with the dispersion of the ID 
antiferromagnetic Heisenberg model 



Uqj ' Ji = 0.6 x — sing 



(9) 



if we include the factor 0.6 as in Eq. (9). The factor is 
independent of J\ for g ^ 1 and at low T. 

These results suggest that at low temperatures there 
is a parameter region where the spin degrees of freedom 
behaves independently from the pseudospin degrees of 
freedom; it is when g > 1 and the gap of the pseudospin 
excitation spectra is large, inside of which there is a spin 
excitation spectra. Thus, we suggest the validity of the 
decoupling of the coupling term of the Hamiltonian as 

L 1 
H ST ^.hJ2 ( T ^+i + H - c -) ( S * ' S *+i - 7) ( 10 ) 



i=i 



with 



(11) 



which leads to the effective Heisenberg-model descrip- 
tion of the spin degrees of freedom of our model. Here, it 
should be noted that in general the factor ^X^T^+H.c.) 
at zero temperature takes the value 1 when g — > and 
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1/2 when g — > oo; thus the above value 0.6 reflects the 
effect of quantum fluctuations of pseudospins, which is 
strong already at g > 1. It should also be noted that 
with increasing temperature the value of the factor de- 
creases to (due to thermal fluctuations), either very 
rapidly when g is small or rather slowly when g is large. 
We find that the essential features of (T^~T~ +1 + H.c.) 
are contained already in a two-site (or dimer) model of 
Eq.(l), a minimum model reflecting the spin-pseudospin 
coupling. This is evident in Fig. 5. 

C. Uniform spin susceptibility 

To see the validity of the effective Heisenberg-model 
description further, in particular for its temperature de- 
pendence, we calculate the temperature dependence of 
the uniform spin susceptibility for the coupled spin- 
pseudospin Hamiltonian. The results are shown in Fig. 6, 
where comparisons are made with the uniform suscep- 
tibility for the system of free spins and with that for 
the ID antiferromagnetic Heisenberg model. We find 
that the temperature k^I '/ Ji at which JiXs{T) shows 
a maximum is lower than that of the ID antiferromag- 
netic Heisenberg model; it becomes lower with decreasing 
the value of g or with increasing the value of J\j Ji. In 
other words, the deviation from the Heisenberg model is 
large when the quantum fluctuation of the pseudospins 
is small, which occurs when g is small or J\ is large. 

Let us analyze the data more precisely. In order to do 
this, we fit the results with the temperature dependence 
of the spin susceptibility of the ID antiferromagnetic 
Heisenberg model, the so-called Bonner-Fisher curve; 28 
i.e., we introduce the T-dependent effective exchange 
coupling constant J e g(T) and determine the values so as 
to fit the calculated uniform spin susceptibility xs(T). If 
the values of J c g thus obtained do not depend on T, it fol- 
lows that the spin degrees freedom of our spin-pseudospin 
model is reduced to a ID Heisenberg model 

L 

Wspin = Jefi J]] Si • Sj+x (12) 

i=l 

at least for the response to the uniform magnetic field. 
The results are shown in Fig. 7. We find that the esti- 
mated value of Jeff(T) is indeed a constant for temper- 
atures below k^T < 0.7 Ji at g = 2. A crossover tem- 
perature T* (= O.7J2) is thereby defined. The effective 
exchange coupling constant thus deduced takes a value 

Jeff ~ 0.6J 2 (13) 

at T < T*; this value is consistent with the value es- 
timated from the dispersion relation of the spin excita- 
tion spectra (see Sec. Ill B). We find that also at g = 4 
the scaling behavior holds up to a higher temperature 
(fceT i$ O.8J2), but with a slightly smaller value of J e ff 
(see Fig. 7 (d)), demonstrating the validity of the effec- 
tive Heisenberg-model description at T < T*. At g = 1, 



however, the temperature region where J c g(T) takes a 
constant value is already very small, although the value 
is still Jog ~ O.6J2 at T ~ K, and at g = 0.5, the value 
of J c ff at T ~ K deviates largely from J e ff = 0.6 J2 (or 
decreases strongly when J1/J2 is large), where the effec- 
tive Heisenberg-model description completely fails. We 
thus find that T* thus deduced is insensitive to J\ , scales 
well with J2, and depends strongly on g (i.e., T* — > 00 
at g — ► 00 and T* ~ at g ~ 1). Note that if we assume 
Eq. (10) the behavior should come from the pseudospin 
fluctuation {T^~T~ Sra + H.c), and actually we find that 
very rough tendency in the parameter and temperature 
dependence is seen in the results for the dimer model (see 
Fig. 5) although the scaling behavior and the presence of 
T* are not seen. 

We note here that the crossover temperature T* 
roughly scales with J2 rather than J\, as seen in Fig. 7. 
One might suppose that it should scale with the size of 
the charge gap: i.e., up to temperatures corresponding to 
the energy of the lowest charge excitations, with which 
the pseudospins can excite, the spin excitations may be 
written in terms of the ID antiferromagnetic Heisenberg 
model. However, as we have discussed in Sec. Ill B, the 
size of the charge gap (if it is defined as a low-energy 
edge of the peak) shows a rather complicated behavior 
and does not simply scale with either J 2 or J\. The 
naive picture thus does not hold. However, since there 
is no other excitations available, the deviation from the 
ID Heisenberg-model description is necessarily due to the 
pseudospin excitations. 



D. Specific heat 

Finally, we present the calculated results for the tem- 
perature dependence of the specific heat C by the finite- 
temperature DMRG methodi^iSi The results are shown 
in Fig. 8, where we find that the curves have a two-peak 
structure: e.g., at g — 2 and J1/J2 = 4, a rather sharp 
peak appears at a low-temperature region k^T ~ O.3J2 
(which scales with J2) and a broad peak structure ap- 
pears at a high-temperature region fee? 1 ^ O.8J2 (which 
scales with Ji), each of which corresponds to the spin 
and pseudospin excitations, respectively. We also find 
that the shape of the low-temperature peak can be fitted 
very well with the temperature dependence of the specific 
heat of the ID antiferromagnetic Heisenberg model with 
the effective exchange coupling constant J c g ~ O.6J2 at 
g = 2; this value is in accord with the value estimated in 
Sec. Ill B. The temperature at which the deviation in the 
fitting occurs is at k^T/ Ji ~ 0.5 — 1 depending on the 
value of g, which is also consistent with the estmate from 
the temperature dependence of the uniform spin suscep- 
tibility. The results thus demonstrate the separation of 
the energy scales between spin and pseudospin degrees of 
freedom and the presence of low-energy 'magnetic' energy 
scale as has been pointed out in Refill 
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IV. SUMMARY AND DISCUSSION 

We have calculated the spin and pseudospin excita- 
tion spectra and the temperature dependence of the uni- 
form spin susceptibility of the coupled spin-pseudospin 
Hamiltonian for the anisotropic Hubbard ladder at quar- 
ter filling by using the QMC method. We have also cal- 
culated the temperature dependence of the specific heat 
of the model by the finite-temperature DMRG method. 
We have first shown that, when the pseudospin quan- 
tum fluctuation is large (g > 1), the dispersion relation 
of the spin exitation spectra of our model at low tem- 
peratures agrees well with that of the ID antiferromag- 
netic Heisenberg model with the renormalized effective 
exchange coupling constant J e g = O.6J2 that is inde- 
pendent of the energy scale of the pseudospin system 
J\. Here, the spin excitation spectra is well inside the 
charge gap, and thus the spin degrees of freedom are sep- 
arated from the charge degrees of freedom. We then have 
shown that the temperature dependence of the uniform 
spin susceptibility of our model is well described again 
by the ID antiferromagnetic Heisenberg model with the 
same effective exchange coupling constant J e g = O.6J2. 
This description is valid up to the crossover tempera- 
ture T* that is related to the pseudospin excitations of 
the system and roughly scales with J 2 unless the quan- 
tum fluctuation of the pseudospins is small (g < 1). We 
have also shown the appearence of the two-peak structure 
in the temperature dependence of the specific heat and 
have confirmed the presence of the low-energy magnetic 
energy scale. We have thus demonstrated the validity 
of the effective Heisenberg-model description of the cou- 
pled spin-pseudospin model for the quarter-filled ladders. 
Then, it follows that the coupling between the spin and 
pseudospin degrees of freedom, which occurs at g < 1, 
leads to the anomalous spin and charge dynamics of the 
system where the spin excitations deviate largely from 
the effective Heisenberg-model description. 

Finally, let us discuss some possible experimental rel- 
evance of our results. The value of the physical param- 
eters appropriate for a'-NaV205 have been estimated in 
Ref£, where we have t\\ ~ 0.14 eV, t± ~ 0.30 eV, and 
Vjl ~ VI ~ 0.8 eV, which lead to J x ~ 1.6 eV, J 2 ~ 0.10 
eV, and g ~ 0.75. We thus find that the real material 
may be in the region of g < 1, where the spin degrees 



of freedom are not completely separated from the charge 
degrees of freedom. The anomalous response of the spin 
degrees of freedom may therefore be expected. We here 
want to point out that the value of J ff estimated from 
the uniform susceptibility observed in experiment (which 
takes the value ~ 600 — 700 K at T ~ K) indeed de- 
creases with increasing temperature^ which is consistent 
with the results of our calculation. Our explanation is 
that the effective spin exchange coupling constant de- 
pends on the fluctuation of pseudospins, the tempera- 
ture dependence of which is strong above the crossover 
temperature (T* ~ K in this material) where the 
spin-charge coupling becomes relevant. The reported^ 
temperature dependence of the nuclear spin-lattice re- 
laxation rate 1/Ti in this material is also interesting in 
the present respect. The measured^ strong temperature 
dependence of the integrated optical spectral weight have 
also been discussed in terms of the destruction of short- 
range spin correlations which occurs as the temperature 
is increasediSi We thus want to point out that observed 
anomalous spin dynamics in the disorder phase of a!- 
NaV2C>5 may possibly be a consequence of this type of 
spin-charge coupling. Because the anomalous charge dy- 
namics has also been noticed in other transition-metal 
oxides^ and some organic systems^, we hope that the 
present study will stimulate further researches on the in- 
triguing interplay between the spin and charge degrees of 
freedom of strongly correlated electron systems with the 
CO instability. 



Acknowledgements 

We would like to thank A. W. Sandvik, T. Mutou, 
N. Shibata, Y. Shibata, and T. Suzuki for useful discus- 
sions on the numerical techniques and T. Ohama for en- 
lightening discussion on the experimental aspects. This 
work was supported in part by Grants-in-Aid for Sci- 
entific Research from the Ministry of Education, Cul- 
ture, Sports, Science, and Technology of Japan. Com- 
putations were carried out at the computer centers of 
the Research Center for Computational Science, Okazaki 
National Research Institutes, and the Institute for Solid 
State Physics, University of Tokyo. 



1 H. Smolinski, C. Gros, W. Weber, U. Peuchert, G. Roth, 
M. Weiden, C. Geibel: Phys. Rev. Lett. 80 (1998) 5164. 

2 H. Seo and H. Fukuyama: J. Phys. Soc. Jpn. 67 (1998) 
2602. 

3 S. Nishimoto and Y. Ohta: J. Phys. Soc. Jpn. 67 (1998) 
2996. 

4 P. Thalmeier and P. Fulde: Europhys. Lett. 44 (1998) 242. 

5 M. V. Mostovoy and D. I. Khomskii: Solid State Commun. 
113 (2000) 159. 

6 M. Isobe and Y. Ueda: J. Phys. Soc. Jpn. 65 (1996) 1178. 



7 T. Ohama, H. Yasuoka, M. Isobe, and Y. Ueda: Phys. Rev. 
B 59 (1999) 3299. 

8 H. Sawa, E. Ninomiya, T. Ohama, H. Nakao, K. Ohwada, 
Y. Murakami, Y. Fujii, Y. Noda, M. Isobe, and Y. Ueda: 
J. Phys. Soc. Jpn. 71 (2002) 385. 

9 D. C. Johnston, R. K. Kremer, M. Troyer, X. Wang, A. 
Kliimper, S. L. Bud'ko, A. F. Panchula, and P. C. Canfield: 
Phys. Rev. B 61 (2000) 9558. 

For a review, see T. Ohama: Bussei Kenkyu (Kyoto) 74 
(2000) 391. 



6 



S. Ravy, J. Jegoudez, and A. Revcolevschi: Phys. Rev. B 
59 (1999) R681. 

H. Nakao, K. Ohwada, N. Takesue, Y. Fujii, M. Isobe, and 

Y. Ueda: Physica B 21-243 (1998) 534. 

A. Damascelli, C. Presura, D. van der Marel, J. Jegoudez, 

and A. Revcolevschi: Phys. Rev. B 61 (2000) 2535. 

C. Presura, D. van der Marel, A. Damascelli, and R. K. 

Kremer: Phys. Rev. B 61 (2000) 15762. 

S. Nishimoto and Y. Ohta: J. Phys. Soc. Jpn. 67 (1998) 

3679. 

S. Nishimoto and Y. Ohta: J. Phys. Soc. Jpn. 67 (1998) 
4010. 

M. V. Mostovoy, D. I. Khomskii, and J. Knoester, : Phys. 
Rev. B 65 (2002) 064412. 

J. Hemberger, M. Lohmann, M. Nicklas, A. Loidl, M. 
Klemm, G. Obermeier, and S. Horn: Europhys. Lett. 42 
(1998) 661. 

A. B. Sushkov, J. L. Musfeldt, X. Wei, S. A. Crooker, J. 
Jegoudez, and A. Revcolevschi: Phys. Rev. B 66 (2002) 
054439. 



20 T. Ohama, M. Isobe, and Y. Ueda: Physica B 329-333 
(2003) 934. 

21 M. Aichhorn, P. Horsch, W. von der Linden, M. Cuoco: 
Phys. Rev. B 65 (2002) 201 101 (R). 

22 M. Cuoco, P. Horsch, and F. Mack: Phys. Rev. B 60 (1999) 
R8438. 

23 D. Sa and C. Gros: Eur. Phys. J. B 18 (2000) 421. 

24 For details, see T. Nakaegawa: Thesis (Chiba University, 
2002). 

25 X. Wang and T. Xiang, Phys. Rev. B 56 (1997) 5061. 

26 N. Shibata, J. Phys. Soc. Jpn. 66 (1997) 2221. 

27 S. Sachdcv: Quantum Phase Transitions (University Press, 
Cambridge, 1999). 

28 J. C. Bonner and M. E. Fisher: Phys. Rev. 135 (1964) 
A640. 

29 R. Amasaki, Y. Shibata, and Y. Ohta: Phys. Rev. B 66 
(2002) 012502. 

30 Y. Shibata, S. Nishimoto, and Y. Ohta: Phys. Rev. B 64 
(2001) 235107. 



7 




aU : aM, u!J 1 



FIG. 2: Dynamical pseudospin structure factor St((7, w) for the coupled spin-pseudospin model calculated at feT = O.IJ2. The 
results at J2/J1 = are for the quantum Ising model. The peak at w = for J2/J1 > in the uppermost panel of (a) and (b) 
is spurious, which is due to the error of the maximum entropy method. 




FIG. 3: Dynamical spin structure factor Ss{q,ui) for the coupled spin-pseudospin model calculated at ksT — O.U2. 
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FIG. 4: Dispersion relations of the spin (open symbols) and 
pseudospin (solid symbols) excitations calculated at UbT = 
O.IJ2. Note that the same data at g = 2 (at g = 1) are plotted 
in (a) and (b) (in (c) and (d)) in different energy scales Ji and 
J2. The dotted line in (a) and (c) is the dispersion relation 
for the quantum Ising model Eq. (5), and that in (b) and (d) 
is the scaled dispersion relation for the ID antiferromagnetic 
Heisenberg model Eq. (9) . 




FIG. 5: Temperature dependence of the pseudospin quantum 
fluctuation (T^T^ +1 + H.c.^ calculated exactly for a dimer 
model of our spin-pseudospin Hamiltonian Eq.(I). 
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FIG. 6: Temperature dependence of the uniform spin suscep- 
tibility xs calculated for the coupled spin-pseudospin Hamil- 
tonian. The solid and dotted curves are the uniform suscepti- 
bility for the system of noninteracting 5=1/2 spins and that 
for the ID antiferromagnetic Heisenberg model, respectively. 




FIG. 7: Effective exchange coupling constant J c s(T) esti- 
mated from the fitting of the calculated uniform spin sus- 
ceptibility to the Bonner-Fisher curvep— Note that the same 
data are plotted as a function of feeT/Ji (left panels) and of 
fcaT / J2 (right panels), whereby a scaling behavior is seen in 
the latter. The arrows indicate the crossover temperature T*. 
The solid lines are the guide to the eye. 
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FIG. 8: Temperature dependence of the specific heat C cal- 
culated for our coupled spin-pseudospin Hamiltonian. 



